LECTURES ON THE ASYMPTOTIC EXPANSION OF A 
HERMITIAN MATRIX INTEGRAL 



MOTOHICO MULASE 



Abstract. In these lectures three different methods of computing the asymp- 
totic expansion of a Hermitian matrix integral is presented. The first one is a 
combinatorial method using Feynman diagrams. This leads us to the gener- 
ating function of the reciprocal of the order of the automorphism group of a 
tiling of a Riemann surface. The second method is based on the classical anal- 
ysis of orthogonal polynomials. A rigorous asymptotic method is established, 
and a special case of the matrix integral is computed in terms of the Riemann 
(^-function. The third method is derived from a formula for the r-function 
solution to the KP equations. This method leads us to a new class of solutions 
of the KP equations that are transcendental, in the sense that they cannot 
be obtained by the celebrated Krichevcr construction and its generalizations 
based on algebraic geometry of vector bundles on Riemann surfaces. In each 
case a mathematically rigorous way of dealing with asymptotic series in an 
infinite number of variables is established. 
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0. Introduction 

The purpose of these lectures is to explain three different methods of calculation 
of the asymptotic expansion of a Hermitian matrix integral. 

The first method is a combinatorial one using the technique of Feynman diagram 
expansion. This method leads us directly to the connection between the matrix 
integrals and the moduli spaces of pointed Riemann surfaces || , |Q , [jllj , |lf| . The 
second method is the classical asymptotic analysis of orthogonal polynomials. It 
allows us to compute the integral explicitly in the special case known as the Penncr 
Model, which is related to the Euler characteristic of the moduli spaces of Riemann 
surfaces. We will see that the values are expressed in terms of the Riemann zeta 
function. Except for this special case, the integral in general reduces to a Selberg 
integral which is not explicitly computable. However, through the fact that the 
Hermitian matrix integral satisfies the KP equations, we give another expression of 
the asymptotic expansion as a t- function of the KP equations. 

The Hermitian matrix integral thus connects three different worlds of mathemat- 
ics: the moduli theory of Riemann surfaces through combinatorics, the Riemann 
zeta function through classical asymptotic analysis, and the theory of integrable 
systems through r- functions of the KP equations. We explain these relations in 
this article, however, no attempt will be made to give any conceptual or geometric 
explanation why the KP equations are related to the topology of moduli spaces of 
pointed Riemann surfaces. 

Riemann's collected work is a great source of imagination to a mathematician. 
The Riemann theta functions were introduced in his monumental paper Theorie 
der Abel'schen Functionen that was published in Crelle's journal in 1857. Two 
years later he published a paper on the prime number distribution where he stud- 
ied the property of the zeta function as a complex analytic function. These papers 
are unrelated, but we note that his proof of the functional equation of the zeta 
function is based on the transformation property of a Jacobi theta function with 
respect to the Jacobi imaginary transform r i— > — 1/r. The Jacobi theta functions 
are the 1-dimensional version of the Riemann theta functions, and the Jacobi imag- 
inary transform is a special case of more general modular transforms in the moduli 
parameters. The coincidental equivalence between the functional equation of the 
Riemann zeta function and the modular invariance of a theta function is myste- 
rious. How much more did Riemann know about the relations between these two 
types of functions? 

In the following sections we explore another relation between these two types 
of functions. The way we will encounter the moduli spaces of Riemann surfaces is 
quite different from Riemann's in the above mentioned paper of 1857. They appear 
very naturally in the asymptotic expansion of Hermitian matrix integrals, which can 
be considered as a kind of generalization of the Riemann theta functions. We know 
that Riemann theta functions associated with Riemann surfaces are characterized 
as finite-dimensional solutions to the system of KP equations jjj , ||, ylj , (l3) . The 
matrix integrals that we will investigate in this article satisfy again the same KP 
equations, though this time they are truly infinite-dimensional solutions 0. 

Using a combinatorial and number-theoretic method, Harer and Zagier Q ob- 
tained a formula for the Euler number of the moduli space of pointed Riemann 
surfaces (defined as an algebraic stack or an orbifold) in terms of the Riemann zeta 
function. Later an analytic method of calculating the asymptotic expansion of a 
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special Hermitian matrix integral was proposed by Penner ^Hj . He discovered that 
the coefficients of the asymptotic series are given in terms of special values of the 
Riemann zeta function. Penner's proposed computation coincides with the formula 
of Harer and Zagier, except for the subtle point of giving an ordering to the set of 
marked points or not. The calculation of the asymptotic expansion of the Penner 
model has been rigorously performed f| . The theorem of Harer and Zagier gives 
an amazing relation between the Riemann zeta function and the Riemann theta 
functions, if we think the latter to be essentially related to the moduli spaces of 
Riemann surfaces. 

We add to this link yet another player: the KP equations. The observation 
|sj that the Hermitian matrix integral is a continuum soliton solution to the KP 
equation is suggestive from the geometric point of view. Soliton solutions repre- 
sent singular Riemann surfaces with rational double points. When we increase the 
number of singularities to continuum infinity the r-function of the soliton solution 
converges to a Hermitian matrix integral that has the information of the Euler 
characteristic of the moduli spaces of pointed Riemann surfaces. We do not know 
why. 

Many explicit formulas for solutions of the KP equations have been established. 
All these solutions are based on the one-to-one correspondence between certain class 
of solutions of the KP equations and a set of geometric data consisting of an arbi- 
trary irreducible algebraic curve, which can be singular as well, and a torsion-free 
sheaf defined on it 0. Let us call a solution to the KP equations transcendental if 
it does not correspond to any algebraic curve. How can we construct a transcen- 
dental solution, then? An answer has been obtained by an accident. It turns out 
that the Hermitian matrix integrals we deal with in this article are transcenden- 
tal solutions of the KP equations. This is closely related to the unexpected sl(2) 
stability condition of the points of the infinite-dimensional Grassmannian of Sato 
12 1 that correspond to the matrix integrals. Again we do not have any satisfac- 
tory explanation why the KP equations, the sl(2) stability condition, and the Euler 
characteristic of the moduli spaces of pointed Riemann surfaces are related. The 
last section is devoted to this topic. 

The organization of the article is as follows. In Section [j] we explain the technique 
of the Feynman diagram expansion through a toy model. A Feynman diagram is 
a kind of graph, but the notion of the automorphism group of a Feynman diagram 
is different from the usual graph theoretic automorphism. This topic is carefully 
treated in this section. Section is devoted to explaining the ribbon graph ex- 
pansion of a Hermitian matrix integral. The mathematical method of dealing with 
asymptotic series in an infinite number of variables is also explained in this section. 
The Penner model is rigorously calculated in Section Ejl following the idea of Q. 
The value we obtain is the Euler characteristic of the moduli spaces of pointed Rie- 
mann surfaces calculated by Harer and Zagier, but we will not go into the moduli 
theory in this article. The third expression of the asymptotic expansion of the Her- 
mitian matrix integral is computed by using the formula for the r-function solution 
to the KP equations in Section ^. This solution is transcendental, which is proved 
in Section ^ from the sl(2) stability condition of the point of the Grassmannian 
that corresponds to the Hermitian matrix integral. The last two sections contain 
our new results, including Theorem |f.2] and Theorem 5.2, which were presented in 
the UIC Workshop in 1997. 
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1. Feynman diagram expansion of a toy model 
Let us start with a simple integral: 

(1.1) / e- x2 / 2 dx = V2^. 



According to Lord Kelvin, a mathematician is one to whom that is as obvious as 
that twice two makes four is to you. However, the usual proof of this formula using 
polar coordinates of a plane is really trivial, and it is hardly a good qualification for a 
mathematician. It is plausible that Lord Kelvin had in mind a proof using functions 
only in one variable and appealing to an infinite product expansion of trigonometric 
functions, that requires reasonably deep knowledge of function theory. 



The integral we consider is a variation of (1.1 



(1.2) Z(t)= I e -* 2 /2 e ^ 4 /4!^ 



/2tt 

We want to know the integral Z(t) as a function of t. Since 



|e t-x /4! | = e Re{tyx /4! 



for every x € R, the integral converges to make Z(t) a holomorphic function in t 
for Re(t) < 0. Unfortunately there is no analytic method to give a simple closed 



formula like (1.1) for (1.2), so we need a different approach. Since a holomorphic 
function defined on a domain is completely determined by its convergent Taylor 
expansion at a point in the domain, we can try to find a convergent power series 
expansion of Z(t). But here again we encounter the same problem, and the only 
thing we can do is restricted to the power series expansion of Z(t) at t = 0. At a 
boundary point of the domain where the function is not holomorphic, there is no 
longer a Taylor expansion, but we still have a useful power series expansion called 
an asymptotic expansion. 

Definition 1.1. Let fi be an open domain of the complex plane C having the origin 
on its boundary, and let h(z) be a holomorphic function defined on f2. A formal 
power series 

oo 

y^^a v z v 

is said to be an asymptotic expansion of h(z) on f2 at z = if 



a v z 



(1-3) I™ * ' XIT =a »'+ 1 



zesi 
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holds for all m > 0. 

If h{z) happens to be holomorphic at z = 0, then the Taylor series expansion of 



h(z) at the origin is by definition an asymptotic expansion. Formula (1.3) shows 
that if h{z) admits an asymptotic expansion, then it is unique. However, we cannot 
recover the original holomorphic function from its asymptotic expansion. Let us 
compute the asymptotic expansion of e 1 / 2 defined on a domain 

(1.4) Q e = {z G C|tt/2 + e < arg(z) < 3tt/2 - e} 

for a small e > 0. Since 

lim —r^- = 

z— o z m+1 

for any m > 0, the asymptotic expansion of e 1 / 2 at the origin is the 0-series. Thus 
the asymptotic expansion does not recognize the difference between e 1 / 2 and the 
0-function. We will use this fact many times in Section |^ when we compute the 
Penner model. This example also shows us that even when h{z) is not holomorphic 
at z = 0, its asymptotic expansion can be a convergent power series. 

To indicate that the asymptotic expansion of a holomorphic function is not equal 
to the original function, we use the following notation: 



A(h(z)) =J2 a vZ v - 



v=0 

If two holomorphic functions h(z) and f(z) defined on f2 have the same asymptotic 
expansion at z = 0, then we write 

h(z)if(z). 

Thus = e 1 / 2 at z = as holomorphic functions defined on the domain Q e . For 
two holomorphic functions f(z) and g(z) defined on Q admitting the asymptotic 
expansions at 0, we have 

A(f(z) + g(z))=A(f(z))+A(g(z)) 

A(f(z)-g(z))=A(f(z))-A(g(z)). 

We note that the asymptotic expansion of a holomorphic function does depend on 
the choice of the domain f2. For example, e 1 ' 2 does not admit any asymptotic 
expansion at z — as a holomorphic function on the right half plane. However, if 

Oi c n 2 , e <9^i n<9ft 2 , 



as in Figure 1.1, and h(z) has an asymptotic expansion on at z = 0, then it also 



admits an asymptotic expansion on ill at z = 0, which is actually the same series. 




Figure 1.1. Domains Oi C ^2 
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We can also define the asymptotic expansion of a real analytic function: if K 
is an open interval of the real axis with as its one of the boundary points and 
h(z) a real analytic function on K, then the same formula (|lj^), replacing fl by K, 
defines the asymptotic expansion of h(z) at z = 0. 

Now let us compute the asymptotic expansion of Z(t) of (1.2) as a holomorphic 
function defined on f2 e = {t £ C|7r/2+e < arg(t) < 3ir/2 — e}. The Taylor expansion 
of the exponential gives 

-^/a t-x*/4i_!*L = r^/ay 1 x 4v f v 

V2¥ J-oo ^(4!)"-«! v^F' 



v=0 

The infinite integral and the infinite sum we have here are not interchangeable. But 
let's just interchange them and see what happens: 

Note that this is a well-defined formal power series in t because the integral 

OO J 

e -x 2 /2 , x 4v . ax 



converges. 



Lemma 1.2. The formal power series (1.5) gives the asymptotic expansion ofZ(t): 



n (4!)°-v! ' V2^ 



E 



Although we cannot get an equality by interchanging the integral and the sum 



because the power series expansion of the integrand of Lemma 1.2 is not uniformly 
convergent on the infinite interval (— oo, oo), at least we obtain a formula which is 
correct in one direction. 



Proof. Using the linearity of the integral, we have 

dx ^ r 



-x z /2 



u=0 v ' 



x 4v ■ t v 



'2tt 



E 

v=0 



(4!)" • u! 



-x 2 /2 _ 4v 



E 

v—m-\- 



i 



=t m+1 



1 (4\) v -v\ 



• a; 411 • t v ■ 



dx 



^ (4!) m + 1 + a ■ (m + l + a)! 

a— s 7 v ' 



dx 



dx 



As long as t stays in £l e = {t E C|7r/2 + e < arg(t) < 3tt/2 — e}, we can divide the 
above expression by t m+1 and take the limit t — > 0, because the integral converges. 
The result is the (m + l)-th coefficient of the asymptotic expansion, which proves 
the claim. □ 
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How can we calculate the coefficient of the expansion? The standard technique 
is the following: 

iv 



e -*72 . x *» 



dx 

7^ 



-x72 . 



4v 



dx 



dx 



I 2 . p*y . 



e -{x- y y/2 _ e ^/2 



/2tt 



7 \ 4f 

dy 



where we have used the translational invariance of the integral (1.1). Note that the 
integration is reduced to a differentiation. All we need now is a Taylor coefficient 
of the exponential function e v / 2 , from which we obtain 



(1.6) 



dx 
7^ 



(Av) 



{2v)\-2 2v 



(4»-l)H, 



where the double factorial is defined by 

(2n - 1)!! = (2n - 1) ■ (2n - 3) • (2n - 5) • • • 5 • 3 • 1. 



The quantity fll.q ) has a combinatorial meaning. Let us denote the differential 
operator d/dy by a dot •. We have 4v dots attacking the fort e y I 2 . Since y is set 
equal to after the operation, if only one dot attacks the fort, the result would be 
just 0: 



A. e y 2 /2 
dy 



y e 



1/73 1 



y=o 



\y=0 



0. 



To obtain a nonzero result, the dots have to attack the fort by pairs: 
2 



,1/72 



y 2 e y /2 | 



oiT/ai 



l B =o 



1. 



y=0 



Noting that the result we get by the paired attack is 1, we conclude that the value 
of the integral (or the differentiation) ([Lq) is equal to 



The number of ways of making 2v pairs out of Av dots 

Av\ fAv ~ 2\ f Av - 4 



(2v)\ 



Av(Av-l) (4v-2)(4u-3) 
" 2 2 
(4u)! 



4-3 2- 1 / . ., 



(2v)\-2 2v ' 

These pairs can be visualized by a diagram like Figure \l.2\ . Let us call such a 
diagram a pairing scheme. Thus ( p_.6| ) gives the number of pairing schemes of Av 



dots. An example of a pairing scheme of 8 = 4 x 2 dots is given in Figure 1.2 
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Figure 1.2. Pairing Scheme 



The coefficient of the asymptotic expansion of Lemma 1.2 has an extra factor 
of 1/(4!)" • vl. How can we interpret it combinatorially? Here enters the idea of 
Feynman diagrams. The 4v dots are grouped into v sets of 4 dots. Let us replace 
each set of 4 dots by a cross, identifying the four dots with the four endpoints of 
the cross. Then the pairing scheme changes into a Feynman diagram, as shown in 
Figure |l.3| , by connecting the endpoints according to the pairing rules. This is an 
example of a graph. We use this word for a CW complex like Figure 1.3 in this 
article. A graph T = (V, E, i) consists of a finite set V of vertices, a finite set E of 
edges, and the incidence relation i of vertices and edges. The number of half-edges 
coming out of a vertex of a graph is called the degree of the vertex. A degree d 
graph is a graph whose vertices have the same degree d. A degree 3 graph is also 
called a trivalent graph. The order of the graph T is the number of vertices \V\ of 

r. 

When we make the Feynman diagram T from a pairing scheme, we consider 
the center of a cross as a vertex and a pairing of dots as an edge of the graph T. 
Figure Ol is thus considered as a degree 4 graph of order 2. 




Figure 1.3. Feynman Diagram 



As a graph we can interchange the v crosses freely, and in each cross we can 
place the four edges in any way we want, as long as the strings are attached. The 
degrees of freedom for these moves are exactly (4!)" • vl. Thus we (tentatively) 
conclude that the v-th coefficient of the asymptotic expansion of the integral Z(t) 
is the number of degree 4 graphs of order v. As an example, let us compute the 
simplest case v = 1. From the above considerations, the number of degree 4 graphs 
with one vertex should be 

(4-1)!! _ 1 
4! ~ 8' 

But this is impossible! What went wrong? 

The number of different pairing schemes of 4 dots is three, as in Figure 1.4. 
When we factored out ( |l.6| ) by (4!)" • vl, we assumed that interchanging the v 
crosses and renumbering each edge of a cross would lead to a different pairing 
scheme that still corresponds to the same graph. In other words, we assumed that 
the group & v x {&i) v acts on the set of all pairing schemes freely, where 6„ denotes 
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the permutation group of n letters, and the product is the semi-direct product of 
two factors with (64)" as its normal subgroup. But as we see clearly from the 
above example, some pairing schemes are stable under the action of non-trivial 
permutations. The isotropy group that stabilizes any of the three pairing schemes 
of Figure |L4] is (Z/2Z) 3 . 



L 



1 



j 



r 



1 1 



Figure 1.4. Pairing Schemes of 4 Dots 



Since our graphs are constructed from pairing schemes, we define the automor- 
phism group in the following manner: 

Definition 1.3. The automorphism group Aut(r) of a graph T is the isotropy 
subgroup Ip of ©i, xi (64)" that preserves the original pairing scheme P. If pairing 
schemes P and P' correspond to T, then the isotropy groups Ip and I pi are con- 
jugate to one another in 6 V x (&4) v - Therefore, as an abstract group, Aut(r) is 
well-defined. 

Remark. Our definition of Aut(r) does not coincide with the traditional graph 
theoretic definition of automorphism. 

The correct interpretation of 1/8 is then l/|Aut(r)|, where T in this case is 
the degree 4 graph with only one vertex, and we denote by |Aut(r)| the order 
of the group. More generally, we can interpret the formal power series ( |l.5| ) as 
a summation over the set of all pairing schemes modulo the group & v x (&4) v , 
which is equivalent to the set of all degree 4 graphs. The contribution of a graph 
r is modified by the weight of l/|Aut(r)|. Summarizing, we have established: The 
asymptotic expansion of the integral Z(t) is given by 

no / 



A 



-x 2 /2 e t-x 4 /4\ 



dx 



E 



y - J - 

^ Aut(r) 

degree 4 graph 



t v . 



r of order 

Since the number of degree 4 graphs with a fixed number of vertices is finite, the 
right hand side of the above formula is a well-defined element of the power series 
ring Q[[t]]. The degree of each vertex of the graph is 4, which is due to the power 
4 in the exponent of the integral. The same argument thus establishes 



Theorem 1.4. The asymptotic formula 



A 



^-x 2 /2t-x 2 >/(2 3 )\ dx 



E 



E 



degree 2j graph 
\ r of order v 



|Aut(r)| 



■t v e 



J 



holds for an arbitrary j > 2. 

We can consider more general graphs with the integral 



(1.7) 



Z{t\,t2, ■ ■ ■ , t2m) — 
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where m > 2 is an integer. The integral converges if t 2m is in the domain Q e of 



(1.4) and determines a holomorphic function on 

(h,t 2 ,--- ,t 2m ) ec 2 " 1 " 1 xn £ . 

We can expand Z(t 1: i 2 , • • • , t 2m ) as a Taylor series in (ti,t 2 , • • • , i2m-i) € C 2m_1 
and as an asymptotic series in t 2m € f2 e at the origin. Fix a value of t 2m £ £ . 
Then 

exp (t 2m -x 2m /(2m)\) 

acts as a uniformizing factor so that the power series expansion of the integrand in 
terms of x converges uniformly on (— oo, oo) for all values of ti, t 2 , • • • , t 2m -\ 6 C. 
Therefore, we can interchange the infinite integral and the infinite sums: 



e i '±x ) ■ ■ ■ exp ( . r*- 1 ) ■ exp < JiH-^ '' ' 



— CO 



C: 



E 



1! J 1 \{2m-iy. J r \(2my. Jy/2n 

^2m ^ ' 2 " , \ __„ 1 +2u 2 H h(2m)ii 2 ^ ^ 



X 



E Wl i.firpi ) "" ( E 



'2m-l 



«i! • (l!)"" 1 / \ ^ n u 2 m-i! • ((2m - 



Ltlo ^ ! • ((2m)!)— J* v^F" 



We already know that 



oc 



— oo 



x 2 /2 [ \^ ^m" 2 '" \ a ,t) 1 +2^2 + --- + (2m)u 2 m ^ X 

v 2m \-({2m)\Y^ ^ 



A e-^ £ . . , 



E 



l2m I e -x 2 /2 x v 1 +2v 2 +--- + {2m)v 2 rr aX 



« 2m !-((2m)!)— V^F 

at t2m = when the top integral is considered to be a holomorphic function in 
t 2m G Q e . Therefore, we have 



(1.8) A 



3 = 

oo 



i y-^ » V2m POO J 

1 l 2m I -x 2 /2 V!+2v 2 -\ h(2m)u 2 n- a2; 



V _ V 

^ ui! • (l!)"i ^ 



e ' x 
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We now apply the Feynman diagram expansion to the above integral. First, we 
have 

J / j \ vi + 2v2-\ h(2m)i>2 m 

c -x 2 /2 x v 1 +2v 2 + --- + (2m)v 2m aX _ [_) e ^ 2 /2 



y=o 

The pairing scheme of the dot diagram has v\ sets of single dot, v-i sets of double 
dots, • • • , and vim sets of 2m dots. Passing to a Feynman diagram, we have a graph 
with Vj vertices of degree j for j = 1, 2, • • • , 2m. Thus 

Theorem 1.5. We have the following asymptotic formula: 



a { f ^ 2/2 ex P [ v = e tA^ ■ n * 

V— Grap t^^ |Ant(T)| 11 

vertices of degree <2m 

where ^(r) denotes the number of vertices of degree j inY . 

Note that the asymptotic series is a well-defined element of the formal power series 
ring 

Q[[*l,*2,--- .*2m]], 

because there are only finitely many graphs for given numbers «i(r), W2(r), • • • , 

f2m(r). 

Definition 1.6. The automorphism group of Y is defined as the isotropy subgroup 
of 

2m 

n e ^ « e ?' 

3=1 

that stabilizes the pairing scheme corresponding to T. 

As before, the definition of Aut(r) as an abstract group does not depend on the 
particular choice of the pairing scheme corresponding to the graph. 

Let us now consider the relation between general graphs and connected graphs. 
Let a v be the number of arbitrary degree j graphs of order v and c v the number 
of con nect ed degree j graphs of order v, where j > 1 is a fixed number. From 
Figure |l.5|, it is obvious that 



(1.9) a v = E 



rii ! • ■ n?.\ • • • ' 

ni+2ri2+3ri3H — =v 

where n,; is the number of connected components with i vertices in a given graph. 
Formula ( |l.9| ) is equivalent to a simple functional relation in terms of generating 
functions: 

oo / oo \ 

(i.io) E a ^ = ex p E c "*i ' 

v=0 \v=l J 

where we use the convention that ao = 1 and cq = 0. 
In a more general case, let v = (v±, V2, • • • , V2m) an d 

2m 

' v II'; 
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Figure 1.5. A Disconnected Graph 

We denote by c v the number of connected graphs with Vj vertices of degree j, where 
1 < j < 2m, and by a v the number of all graphs with Vj vertices of degree j. Then 
we have 

a v t v = exp Cv ^ V ^ 

where the n-th term of the right hand side counts the number of graphs consisting 
of n connected components. The factor l/(n\) means that we can permute the n 
connected components without changing the original graph. 

The case we are considering is slightly more complicated because the generat- 



(1.11) 



ing function we have in Theorem 1.5 counts the number of graphs with weight 
1 /|Aut(r)|. The automorphism group of a graph T consisting of n connected com- 
ponents T±, ■ ■ ■ , T n is the semi-direct product 

n 

(i.i2) Aut(r) = s„ m TT Aut(r\,-). 



Therefore, we have 



1 n 
7T\ IT 



|Aut(r)| n\ 11 | Author 
Note that the right hand side is the product of n factors following the key factor 



l/(n!). It shows, therefore, that the exponential formula ( 1.11 ) connecting con- 
nected graphs and general graphs also holds in the case we are investigating. Thus 
we have established 

Theorem 1.7. The asymptotic series involving only connected graphs is given by 
log A 




1 



2 m 



IK 



^ lAutrni 11 1 

Connected graph T with j— 1 

vertices of degree <2m 



We note here that the function log is not considered as an analytic function. It is 



applied to the formal power series appearing in the right hand side of Theorem 1.5 
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as the inverse power series of exp(z) denned by 

oo 

iog(i-*) = -£- 



z 

n 

n=l 



2. Matrix integrals and ribbon graphs/fatgraphs 

Let Tin denote the space of all n x n Hcrmitian matrices. It is an rt 2 -dimensional 
Euclidean space with a metric 




v/trace(X - Y") 2 , X, Y e H n . 

The standard volume form on 7i„, which is compatible with the above metric, is 
given by 

d(j,(X) = dxn A dx22 A • • ■ A dx nn A I f \ d(ReXij) A d(Imxi 

\i<j 

for X — [xij] € Ti. n . It is important to note that the metric and the volume form 
of Tin are invariant under the conjugation X i — ► UXU^ 1 by a unitary matrix 
U G U(n). The main subject of this section is the Hermitian matrix integral 

(2.1) Z n (t,m) = J exp ^-itrace(X 2 )^ exp 
where 

(2.2) N = J exp ^-itrace(X 2 )^ dfi(X) = 2"/ 2 • tt™ 2 / 2 

is a normalization constant to make Z n (0, m) — 1. We note that Z n (t,m) is a 
holomorphic function in (£3, t±, ■ ■ ■ , t2m-i) G C 2m_3 and 

t 2 m en e = {te C|tt/2 + e < arg(t) < 3tt/2 - e} 

(e > 0), because the dominating term trace(X 2m ) is positive definite on TL n . Thus 
we can expand Z n (t,m) as a convergent power series in t$, £4, • • • , 42m-i about 0, 
and as an asymptotic series in t2 m at = 0. 

We also note here that we do not include the t\ and t 2 terms in the integral 
because of our interests in topology, which will become clearer as we proceed. From 
the point of view of graphs, we do not allow degree 1 and 2 vertices in this sec tion . 



Corresponding to the fact that the integral ( 2A ) has richer structure than 



1. 



the Feynman diagrams appearing in the asymptotic expansion of Z n (t,m) have 



more information than just a graph as in Theorem 1.5. As we are going to see 
below, the new information we have from the Hermitian matrix integral is that the 
graph is drawn on a compact oriented surface. Suppose we have such a graph drawn 
on an oriented surface, as in Figure [2~l| . 

Locally at each vertex of the graph, the orientation of the surface gives rise to a 
cyclic order of the edges coming out of the vertex, as shown in Figure G. 



A graph drawn on a surface thus gives a graph with a cyclic order of edges at each 



vertex. An example, that is corresponding to Figure 2.1, is shown in Figure 2.3 



Note that two circles are reversed in Figure 2.3, corresponding to the fact that 



two edges of the graph of Figure 2.1 go around the back side of the surface. 
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Figure 2.1. Graph on a Surface 




Figure 2.2. Cyclic Order of Edges 




FIGURE 2.3. Ribbon Graph 



Conversely, suppose we have a connected graph with a cyclic order of edges 
assigned to each vertex. To indicate that we have the extra information of cyclic 
order at each vertex, we use T r n> and distinguish it from the underlying graph T. 
We can construct a compact oriented surface C(T r n,) canonically such that the 
graph r is drawn on it, as follows. First, the graph around each vertex can be 
drawn on a positively oriented plane that is compatible with the cyclic order. Next 
we fatten the local part of the graph into a crossroad of multiple intersection. The 
orientation of the plane defines an orientation on each sidewalk of the crossroad, as 
in Figure 2.4 



OIIO 



olto 



Figure 2.4. Oriented Crossroad 
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The roads are connected to the other parts of the graph, with matching ori- 
entation on the sidewalks. Then we obtain an oriented surface with boundary. 
Figure |2.5| shows such a surface with boundary. 






Figure 2.5. Ribbon Graph and Surface with Boundary 



Let b(T r ib) denote the number of boundary components of this oriented surface 
(= fattened graph) made out of T r n,. From the construction, each boundary com- 
ponent has a unique orientation compatible with that of the fattened graph. Thus a 
boundary component is indeed an oriented circle, which we also call a boundary cir- 
cuit. So we can attach an oriented 2-dimensional disk to each boundary component 
of the fattened graph to construct a compact oriented surface C{Y r u,)- 

Definition 2.1. A ribbon graph (or a fat graph) is a graph with a cyclic order of 
edges assigned to each vertex. 

The ribbon graph r r jb of Figure has only one boundary component, and the 
resulting compac t su rface C{T r ib) is the 2-torus on which the underlying graph T 
is drawn (Figure 2.6). 



Figure 2.6. 3-valent Graph on a Torus 

We have shown that every graph drawn on an oriented surface is a ribbon graph, 
and conversely, that every connected ribbon graph T r ib gives rise to a canonical 
compact oriented surface C{Y r ib) on which the underlying graph is drawn. The 
attached boundary disks and the underlying graph V give a cell-decomposition of 

c(r ri b). 

Lemma 2.2. Let T r n, be a connected ribbon graph with vertices of degree > 3, and 
Vj(T) denote the number of vertices of the underlying graph T of degree j. Then 
the genus g(C(T r it>)) of the canonical oriented surface C(T r n,) associated with T r ib 
is computed by the following formula: 

2 - 2g(C(T rlb )) = ]T Vj(T) - \ £ j ■ v^T) + s(T rib ), 

j>3 j>3 
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Proof. The total number of vertices of the cell-decomposition is given by X)j>3 v j(^)- 
Since each edge is bounded by two vertices (possibly the same) , the number of edges 
is given by ^ X^>3 3 ' v j CO- By construction, b(T ri i,) is the number of 2-cells. Thus 
the Euler characteristic of a compact surface gives the above formula. □ 

To see how ribbon graphs appear in the matrix integral, let us consider a simple 
example: 



1 



exp trace(X ) • exp -traceX 



t 



N 



Using the same argument as in Lemma 1.2 we can prove the asymptotic formula 

dfj,(X)\ 



A exp ^— itrace(X 2 )^ • exp ^-traceX' 



N J 



E 4Td j H 6XP (-^ trace (^ 2 )) ' (traced 4 ) 



N 



We need another matrix Y — [yij] £ H n and a differential operator 



d_ 



d 



to compute the asymptotic expansion of the integral. 
Lemma 2.3. For every j > and v > 0, we have 



(2.3) 



trace 



— V ] e tra " 1 v 1 
dY 



= (traced)" . 



Y=0 



Proof. Suppose that Y and X are both arbitrary complex matrices of size n. Then 
for each j > 0, we have 



trace 



d_ 

dY 



tracc(X t -F) 



E 



Y=0 

8 d 



=1 dy ili2 dy i2l3 dyi 3il 

71 

E 



exp 



xm ■ yu 



Y=0 

= tr&ceX 1 . 



Repeating it v times, we obtain the desired formula (2.3) for general complex ma- 
trices. Certainly, the formula holds after changing coordinates: 



(2.4) 



-ltVi 



for i < j 
for i < j 



where Uij and Wij are complex variables. Since ( |2.3| ) is an algebraic formula, it 
holds for an arbitrary field of characteristic 0. In particular, (2.3) holds for real Uij 



and Wij, which proves the lemma. 



□ 
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Therefore, we have 
exp 



hr&ce(X 2 )^j ■ (traced 4 )" 



dfx(X) 
N 



exp ^~itrace(X 2 )^ • ^ 



Y=0 



dfi{X) 
N 



trace 



trace 



E 



dY 



d 



exp ( — -trace(A — Y 



Y=0 



dn(X) 
N 



3 l/2traceY' 



dY , 

d d d d 
dyij dy jk dyu dy ei 



exp 



yijVji 



Y=0 



The only nontrivial contribution of the differentiation comes from paired deriva- 
tives: 



_d d_ 

dy^ dyke 



exp 



1,3 



ViiVji 



d 

dy i: 



-yik 



lit ■ Ojk- 



Y=0 



If we denote by ».y the differential operator , then we have a pairing scheme of 
4w dots as before, and the pairing of two dots and »ki contributes Sn ■ Sjk- Thus 



trace 



(2.5) 



_d_ 

dY 



E 



l/2traceY' 



Y=0 



E E 



2i ,j\ ,fci,^i=l i v ,j v ,k v ,£ v —l All pairings 

P of 4u dots 



n 



V ■ Ojfc 



All paired dots 



A symbolic description of the contribution of parings is given in Figures 2.7. 



r~i r~i 



jk kl 



SikSjjSkiSee 



jk 



kl 



— SuSjkSjiSki 



jk 



kl 



I_J 



— SuSjeSjgSkk 



Figure 2.7. Pairing Contribution 



An interpretation of Figure 2.7 in terms of Feynman Diagrams was introduced 
by 'tHo oft [m| . The set of four indexed dots »ij »jk *kt*U is replaced by a crossroad 
(Figure 2 



Since 



d 



is different from •,• 



the different roles of the indices are 



represented by an arrow. If is connected to »jk, then it gives a contribution of 
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I k 



Figure 2.8. Indexed Crossroad 



5ik • Sjj. 'tHooft visualized this situation graphically by making a crossroad loop 



(Figure 2.9) 




Figure 2.9. Crossroad Loop 



Note that the orientation of the sidewalks of this crossroad loop is consistent. 
Thus we obtain a ribbon graph, as we expected. 

The passage from the pairing scheme to a ribbon graph has again some redun- 
dancy. In Section [j], the permutation group &j appeared for a vertex of degree j. 
This is due to the fact that a scalar monomial X1X2 • ■ ■ Xj is invariant under the &j- 
action. In the case of matrix integrals, a monomial is of type trace(XiX2 • • • Xj), 
which is invariant under the action of the cyclic group Z/j'Z, but not under the full 
symmetric group &j . This is the origin of the appearance of the extra cyclic order 
of the edges at each vertex. 

Definition 2.4. Let P be a pairing scheme of indexed dots and T r ib the corre- 
sponding ribbon graph. Then the group 

n (©^(D x (wr (r) ) 

acts on the set of all pairing schemes. As before, we define the automorphism group 
of a ribbon graph T r u, to be the isotropy subgroup of the above group that fixes 
P. As an abstract group, Aut(IYi6) does not depend on the choice of the pairing 
scheme of indexed dots corresponding to T ri i,. 

One more difference between the matrix integral and the integrals considered in 
Section |l| is the appearance of the size of matrix in the calculation. To illustrate 
this effect, let us continue our consideration of the degree 4 case with one vertex: 
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4-1! 



3 -l/2trace(X 2 ) trace(x 4 ) ^) 



(2.6) 



1 " 

7 E (SikdjjSkiSu + SuSjkSjiSu + SuSjiSjiSkk) 



-in 3 + n + n 3 ) = -n 3 + -n. 
4 V 7 2 4 



As shown in Figure 2.10 there are two degree 4 ribbon graphs of order one. The 
one on the left has the automorphism group Z/2Z, while the second has Z/4Z. 





FIGURE 2.10. Degree 4 Ribbon Graphs with 1 Vertex 



We also note that the exponent of n in (2.6) is exactly the number of boundary 
components of the ribbon graph which is considered as a surface with boundary. 
For every v > 1, we now have 



1 



4" • v\ 



-l/2t r ace(X=)( trace(x 4 ) ^^l 



1 



^ |Aut(r ri 6*)| 

degree 4 ribbon graph 
I" Vib of order v 



The same argument that we used to prove Theorem 1.5 works and we have 



Theorem 2.5. The asymptotic expansion of the Hermitian matrix integral (2.1) 
is given by 



A 



J exp ^-itrace(X 2 )^j exp ^trace ^ — X J ^ 



dn(X ) 
N 



E 



MT„ b ) 



Ribbon graph V r ib with 
vertices of degree 3,4,-- - ,2m 



|Aut(r 



rib ) 



2m 

■n^ (r) 

3=3 



where b(T r ib) denotes the number of boundary components of the ribbon graph T r ib, 
and Vj(T) the number of degree j vertices in the underlying graph T. 

Here we note that for given values of ^(r),--- ,W2m(r), the number of ribbon 
graphs is finite. Thus the above asymptotic series belongs to (Q[n]j [[t 3 , t±, ■ ■ ■ , t 2m \\ ■ 
The relation between connected ribbon graphs and arbitrary ribbon graphs are the 
same as in Section |l|. In particular, since (1.12) also holds for ribbon graphs, ap- 
plication of the logarithm gives us 
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Theorem 2.6. 



log A I J exp ^-itrace(Jf 2 )^ exp | trace ^ — X j 



3=3 ) j 

. 2m 

4? |Aut(r„ b )| 11 tj 

Connected ribbon graph r r a> J = 3 
with maximum degree 1m 

This formula is particularly useful, because we are interested in connected Riemann 
surfaces and only connected ribbon graphs give rise to connected surfaces. Using 
Lemma 2/2, we can rearrange the summation in terms of the genus of a compact 
oriented surface and the number of marked points on it: 

log A (Z n (t,m)) 

( \ 



(2.7, __ 



3>0,s>0 
2-2g-s<0 



„ 2m 
\' " TT +»i ( r ) 



ik 



^ |Aut(r„- b )| 11 j 

Connected ribbon graph I~V^ J=3 
with vertices of degree 3,4,- ■■ ,2m, 

\ X(r)=2-2 9 - s ,6(r rii ,)=s / 



where X(r) denotes the Euler characteristic of the underlying graph F. Note that 
Lemma 2.2 implies that 

2 - 2g(C(T nb )) - b(T rib ) < 

for every ribbon graph T r n,. Let v(T) and e(T) be the total number of vertices and 
edges of the graph L, respectively. Then 

fx(r) = v (r)-e(r) 

(2.8) I «(T) = w 3 (r) + v 4 (T) + ■■■ + v 2m (T) 

(e(T) = |(3-» 3 (r)+4.«4(r) + .-. + 2m-w ain (r)), 

because the vertices of T have degree in between 3 and 2m. Thus for every fix ed a 
and s, the second summation of (2/?) is a finite sum, which again shows that ( [2.7|) 
is an element of the formal power series ring 

(Q[n])[[t 3 ,U,--- Mm]]- 

The number g is of course the genus of C(r rio ). The topological type of the rib- 
bon graph T ri b is the same as the compact surface C(r ri b) minus b(r rio ) points. 
The number of boundary components becomes the number of marked points of a 
Riemann surface in later sections. 

Let (Q[n]) [[t 3 , t 4 , • • • ]] be the formal power series ring in infinitely many variables. 
The adic topology of this ring is given by the degree 

degt,*=j, j>3 

and the ideal 3j(t) of (Q[n]) [[t 3 , ti, ■ ■ ■]] generated by polynomials in t 3 ,t^, ■ • • of 
degree greater than j, with coefficients in Q[n]. We have a natural projection 

TTj : (Q[n])[[t 3 ,U,---}] — ► (Q[n])[[t 3 ,*4,---]]/3fj(*) = (QW)[[*3.--- MftM 
For each fixed j, the projection image 

7T d (log A(Z n (t,m))) e (Q[»])[[«a,t4,--]]/3j(*) = (Q[n])[[* 3) • • • ,*j]]A-(*) 
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is stable for all 2m > j. Since 
(QN)[[i a) *4,' 



\im(Q[n])[[t 3 ,t 4 ,---}}/3 J (t) 



and 



{7r 2m (log^(Z„(t,m)))} ? 



defines an element of the projective system, it gives a well-defined formal power 
series in infinitely many variables. We denote it symbolically by 



(2.9) 



lim \ogA(Z n (t, to)) 



= {TT 2m (\ogA(Z n (t,m)))} m > 2 e (Q[n])[[t 3 ,U, ■■■]}. 
Going back to the Feynman diagram expansion (pj), we have an equality 



(2.10) 



lim \ogA(Z n (t,m)) 

n. — >oo 

= E 



9>0,s>0 
2-2g-s<0 



E 



|Aut(r„ 6 )| 11 3 

Connected ribbon graph T r ib 
with vertices of degree >3, 

\ X(r)=2-2 S - S ,6(r ri6 )=s / 



as an element of (Q[^])[[<3, *4, • ■ •]]■ For each fixed g and s, the maximum possible 
valency of the ribbon graphs in the second summation is Ag + 2s — 2. To see this, 
let r be a graph with the largest possible degree £. Since the Euler characteristic 
of r is given by 2 — 2g — s = v(T) — e(T), the degree becomes maximum when T 
has only one vertex. Thus 

2 - 2g - s = 1 - -I. 



This shows us that the right hand side of (2.10) does not have any infinite products. 



3. Asymptotic analysis of the Penner model 

There are no known analytic methods to compute the matrix integral Z n (t, m) 
for general m. It is therefore an amazing observation of Penner that at the limit of 
to — > oo a certain specialization of Z n (t, to) is actually computable. In this section 
we study the Penner model and calculate its asymptotic expansion analytically. 

The specialization Penner considered is the substitution 

(3.1) *i = -(Vi) i_2 , i = 3,4,---,2m 

in the matrix integral Z n (t, to), where \fz is defined for Re(z) > 0. The condition 

tt/2 + e < arg(i 2m ) < 3tt/2 - e 
for t-im translates into the condition 

(3.2) | MgW |<_JL_. 
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Thus we have a holomorphic function 

P n (z,m) 



(3.3) 



J H 6XP ( _ ^ trace ^ 2 )) exp ^ _ E 



exp 



(v^) 



i-2 



2m 



J2 ^i^trace(^) 



-trace(X J ) 



d/j.(X) 
N 



i=2 



N 



defined on the region of the complex plane given by (3.2). 




Figure 3.1. Wedge-shape Domain 



We note that the domain (3.2) still makes sense as the positive real axis when 
we take the limit m — > oo. The quantity N is the same normalization constant as 



(2.2) 



The asymptotic expansion of (3.3) at z = can be calculated by making the 
same substitution ( |3.1| ) in Theorem [2.5|. Taking the logarithm, we obtain 

(3.4) 

log A (P n (z,m)) 



E 



3>0,s>0 
2-2g-s<0 



E 



Connected ribbon graph T r ib 
with vertices of degree 3,4.--- ,2m, 
\ X(r)=2-2g- S ,6(r rii) )= S 



(-i) e ( r ) 
|Aut(r 



rib I 



(-*) 



2g+.s-2 



where we used ( ^.8[ ) to compute 

2»i 



7=3 



= (-1) 



«(r)-, e (r)-«(r) 



(-l)^(-z) 



-X{T) 



Note that the right hand side of ([TJ) is a well-defined element of For 
every v > 0, the terms in log A (P n (z,m)) of degree less than or equal to v with 
respect to z are stable for all m > v + 1 . Again by the same argument we used in 



ASYMPTOTIC EXPANSION OF MATRIX INTEGRAL 



2li 



Section || we can define an element 

lim log A(P n (z,m)) e (QM) [[*]]. 

m — »oo 

Thus we have an equality 
(3.5) 

■g^traceCX^^ 
/ 



lim log .4 



exp 



N 



E 



g>0,s>0 
2-2g-s<0 



E 



Connected ribbon graph T r 
with vertices of degree >3. 
\ X(r)=2-2g-s,6(r riil )= S 



(-i) e < r ) 
|Aut(r ri6 )| 



(-*) 



2g+s-2 



J 



as a well-defiiied element of (Q[n]) [[z]]. We recall that in ( 2.1C ) we proved that the 
number of ribbon graphs in the second summation for fixed g and s is finite. 

Let us now compute lim™—,,^ log .4 (P n (z, m)). The standard analytic technique 
to compute the Hermitian matrix integrals is the following formula. Let f(X) be a 
function on X G Ti n which is invariant under the conjugation by a unitary matrix 
U € U{n): 

f(X) = flu- 1 -X-U) = f(k , fei,---, fc„_0, 

where fco, fci, • • • , fc n -i are the eigenvalues of the Hermitian matrix X. If f(X) is 
integrable on Ti. n with respect to the measure dfi(X), then 



(3.6) 
where 
(3.7) 
and 



H„ 



f{X)d»(X) = c(n) • / /(fco, h, • • • , fc„_ 1 )A(fc) 2 dfc * 1 • • • dk n . 



c(n) 



r n(n-l)/2 



n\ ■ (n -!)!••• 2! •!!' 



A(fc) = A(fc ,fci — =det 





A 


fc 


fc 2 






1 


fci 




. . k 1 ^ 1 


det 


1 


fc 2 


fc 2 


un—l 
■ ■ ft 2 




Vi 


fc«— i 


fc 2 





is the Vandermonde determinant. The proof of ( p.6[ ) goes as follows: 

o 

Let 7i„ denote the open dense subset of Ti n consisting of non-singular Hermitian 
matrices of size n with n distinct eigenvalues. If f(X) is a regular integrable function 
on TL n , then 



f(X)dfi(X) 



H„ 



f(X)d»(X). 



We denote by R™ the space of real diagonal matrices of all distinct, non-zero eigen- 

o 

values. Here again integration over R" is equal to integration over R™. Since every 
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Hermitian matrix is diagonalizable by a unitary matrix, we have a surjective map 



U(n) xl"9 U, 



fc 



u 



k 



kn—l 



■ u- 1 e H r 



The fiber of this map is the set of all unitary matrices that are commutative with 
a generic real diagonal matrix, which can be identified with the product of two 
subgroups 

T n • W n C U(n), 

where T n C U(n) is the maximal torus of U(n), and W n C U(n) the group of 
permutation matrices of size n. Note that 

dim U (n) = dim H n = n 2 , dim T n = n. 

Therefore, the induced map 

h : U(n)/T n x W 1 — ► H n 

is a covering map of degree \W n \ = n\. We need the Jacobian determinant of h. 
Put 

~k 



X = [ Xij ] = U ■ 



U- 1 &H n , 



and denote 
Then 

dX = dU ■ 
-U ■ 



h) 



k n -i 



kn—l 



dX = [dxij~\ . 

dko 

■U^ + U- ' 
■U- 1 -dU-U- 1 



dkn-l 



u- 1 



Ti- 



ll- 



= U ■ 



dk 



dk, 



dko 



dko 



dk n -! 



II- 1 + [dU ■ U-\X] 
ko 

IJ- 1 ■ dU, 



k n — l 



■ u- 1 



+ [(kj - fei)dw y ] • U-\ 



where 



rffc„_i_ 

If' 1 -dU= [duji 3 ] 
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which is a skew Hermitian matrix. In terms of the above expression, we compute 
dfj,(X) = dko A • • • A dfc„_i A f\(kj - h) 2 'Re(du>ij) A $s(duJij) 



\i<j 



Thus the integration on Tt n is separated to integration on U (n) jT n and K" . Let 
c(n) = — A Re(dujij) A 3?(c£wy). 

Then we obtain 

f(X)d(i(X) = c(n) i A(k) 2 f(k , • • • , fc„-i)dfco ■ ■ • *n-i- 

For computation of c(n), we refer to, for example, Bessis-Itzykson-Zubcr 
Using formula (|3.6[), we can reduce our integral to 



PnM= c Mf A ( ^n exP -£^ J '" 2 



»=0 \ \ j=2 J 



At this stage, one might want to compute 



2m 



lim exp I — (Vf) u J _ ex p f i log(l — -Jzh N 



= (1 - \Azfc,) 17 " ' e fcl/V5 
= z 1 / 2 -e 1/z -x 1/2 -e- :E , 

where 

1 - y/zh 
x = . 

z 

Since the above function in x is proportional to the Laguerre potential, one might 
expect that the integral becomes computable. However, such a substitution requires 
a very careful treatment. First of all, we have to justify the limit m — > oo taken 
inside the integral over the whole space. Secondly, the integral with respect to ki 
is for the entire real axis, which translates to an integral in x again on the entire 
real line. Since the Laguerre potential is not integrable for negative x, the above 
formal computation cannot be justifiable inside the integral sign. What should wc 
do, then? 

The following is our key idea to compute the Penner model. 

Theorem 3.1 (0)- Let J v (z) — z v ■ C[[z]] denote the ideal of C[[z]} generated by 
z v , and 

tt v : C[[z}\ — > C[[z]]/3 v (z) 

the natural projection. For an arbitrary polynomial p{k) € C[fc], consider the fol- 
lowing two asymptotic series: 

a{z, m) = A\ I p{k) • exp I - V ^? ti I dk I e C[[z}] 



poo I 27n 

/ P( k ) ' exP - X/ 

/- °° V j=2 
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b{z) = ALrz{ez) l ' z J™ p(- 



x 



as z — > +0 with |arg(z)| < 2 m-2 > and 

1 — zx 
Hz 

as z — > +0 with z > 0. Then for every m > 2, we have 

ir m (a{z,m)) = n m (b(z)) 
as an element of < C[[z]]/3 m (z) . In other words, 

2m , /—\j-2 



1/z ■ e~ x ■ dx ] e C\\z 



lim A I I p(k) • cxp ( - V ^ J fcJ ] dfc 
— \ U J J J 

= .4 (Vm 1 ^ jf °° P (^0) ■ xl ' Z ■ e ' X ■ dx ) 
holds with respect to the 3 m (z)-adic topology ofC[[z\\. 

Remark. The above integrals are never equal as holomorphic functions in z. The 
limit m — > oo makes sense only for real positive z, and the equality holds only 
asymptotically. 

Proof. Putting y = \fzk, we have 



/oo I 2m 

p(k) ■ exp - ^ 
V J=2 



k 3 dk 



J 



■ exp 



dv{y,m), 



where 

m) = -L . „ . exp M E J j ^ 

Let us decompose the integral into three pieces: 

/oo y> — 1 />1 />oo 

dv{y 1 m) = I dv{y,m)+ I du(y,m) + / di/(y,m). 
-oo J ~oo J — 1 «/ 1 

Note that the polynomial Y^j=2 ~f °f degree 2m takes positive values on the in- 
tervals (— oo, — 1] and [l,oo). Since p(k) is a polynomial, it is obvious that the 
asymptotic expansion of the first and the third integrals of the right hand side of 
(3.S) for z — > +0 with z > is the 0-series. Therefore, we have 

A f 1 

dv(y,m) = / dv(y,m). 
i J -i 

On the interval [—1, 1], if we fix a z such that Re(z) > 0, then the convergence 

2m . 



lim exp --E^ -(i-y) 17 



z * — ' 7 

3=2 J 
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is absolute and uniform with respect to y. Thus, for a new variable t = 1 — y, we 
have 



m — >oo 



lira / dv(y, m) 



' / p\ —f^r \ % e 



= J- e V* f p ( X -=l ) t^e-^dt - ^e 1 '* [ p ( ) tL/'e-'t'dt. 
Vz Jo \ Vz J V* h \ V* J 

This last integral is 

Since 1 + logt — t < for t > 2, the asymptotic expansion of this integral as 
z — > +0 with z > is the 0-series. Therefore, since the integrals do not depend on 
the integration variables, we have 

Mm .4 I / »(/,-)■ <-xi» I - > ' ^ J I ,//,• 



m — >oo 



(/■oo / 2 m 

y p(fc)-expf-^ 



J 



o 



'2 



as a formal power series in z. This completes the proof of Theorem. □ 



By applying Theorem 3.1 for each fcj, we obtain 



(„ n— 1 / 2m 

/ A(fc) 2 • TT exp - V 
®" to V ^ 



(3-9) =A\Ui£W*) J -J A^i^j • If,-, 1 •« "</,-, 



i=0 

/•oo />oc n— 1 

= X( (v^V^Y***^ / •••/ A(xf • TT x\ ,z e- Xi dxi 
V Jo Jo l=a ) 

where we used the multilinear property of the Vandermonde determinant. We 
can use the standard technique of orthogonal polynomials to compute the above 
integral. Let pj (x) be a monic orthogonal polynomial in x of degree j with respect 
to the measure 

d\(x) = x 1/z e-' x dx 
defined on K = [0, oo) for a positive z > 0: 



K 



Pi{x)pj{x)d\{x) = 5ij || pj(x) || 2 
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Because of the multilinearity of the determinant, we have once again 

A(x) = det (x{\ = det (pj(xi)) . 

Therefore, 



K" 



A(x) 2 d\(xo) ■ ■ ■ d\(x n -i) 

det (pj(xi)) det (pj(x i )dX(x l )) 



K" 



(3.10) 



/n — 1 n — 1 

^2 X! si S n ( (J ) si g n (' r ) II^M^) W_Pr(i){Xi)dX{xi) 

tr£6„ t€6„ i=0 i=0 

n—1 « 

= X! X! sign(<r)sign(r) JJ / p CT(i) (a;)p T(i) (x)dA(x) 
ffes„ re6„ 1=0^ k 

n — 1 

= X! si S n ( CT ) 2 II / P ( r(i)( a: )Pa( l )(a;)rfA(a:) 



o-ee„ 

n-l 

n\ 1] II P*( 



For a real number z > 0, the Laguerre polynomial 



m+lM (-l)J' . (-1)' 



m-j J 3 



of degree m satisfies the orthogonality condition 



(3.11) 

Thus we can use 
(3.12) 



L\' z {x)Ly z {x)e- x x 1 ' z dx = 5 ij U+ ^ /z)l 



i {x) = {-l) i -i\-L\' z {x) 



for the computation. From (3.9)— ( 3.12), we have 



(„ n—1 / 2m 

»=0 \ J=2 



i-2 



-ft? dfe 



(3.13) 



1 



■ „4 (ez) " • 2 "2 • n 
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Applying (3.6) and ( 3.13 ) to (3.5), we conclude 



lim log .A ( 4 / ex P I ~E 

' n ^°° \ N Jn n \ j~2 



J 

n-1 



trace(X J ) dfi(X) 



log A 



1 „(„-!) „ „2 f 1 

7T^— . ( ez )- • • [[ I -1 + - 
i=0 ^ 



iV 



1 



(ez)" • z~ • T - 



2 

n n— 1 



(3.14) 



z^ 

n n, n. 2 , , . / / 1 

const H 1 — log z H log z + n log .4 F 

z z 2 V \ z 



n 

i=0 



i + - 

z 



^(n-i) log 



1 + 2Z 



i=0 



const H 1 — log z — — log z + n log ^4 I T ( — 



n n 



z z 



E 



(-l) r 



5> 

vi=0 



Let us recall Stirling's formula: 



(3.15) log^frfi)) =-ho g z-± + ±logz + jr 



where b r is the Bernoulli number defined by 



^ 2r(2r- 1) 

r— 1 v 7 



2r z2r -! 



const, 



_ = V -x r 
1 ^r! 



r=0 



We are not interested in the constant term (the term independent of z) of (3.15) 
because the asymptotic series in question, ( |3.5| ), has no constant term. We can see 
that substitution of ( [3.15 ) in (3.14) eliminates all the logarithmic terms as desired: 



lim log.4 — / exp ( - V 



(Vi) 



3-2 



-trace(X : 



)j dn(X)j 



E 



>2r 



v2r- 



2r(2r- 1) 



" + E 



r-l 



^(n-i)i 



r=l 



Let 



<f) r (x) = 

9=0 

denote the Bernoulli polynomial. Then we have 

4>r+i(n) 
r+1 



E^ = 



Thus for r > 0, 



:!() 
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E, -N-r n 4>r+l{n) 4>r+2{n) 
[n — i)i = 
V ; r + 1 r + 2 

4=0 



1 /V + lV ^,_„ ^ 1 /r + 2 



y^_ r^u . n ^-«_y ( 

rl(l-g) 



6„ • n r+2 -« 



/ . / „ \i u q n <Jr+l '*• 

^q\(r + 2- q)\ 



Therefore, we have 



(3.16) 



br.n.^-i + y til 

- - E b» • - ■ ^ + EEnr^H^ ■ - r+2 ^ ■ « r 

i 

■ 7 

(r + l)(r + 2) 

(2 g )!(r + 2-2 g )! /l " 



r=l v ; r=l \i=0 / 



CO 1 OC - 

2r ^ K ' r(r + l) 



r=2 ?=1 



It is time to switch the summation indices r and q to g and s as in (3.5). The 
first sum of the third line of (3.16) is the case when we specify a single point on 
a Riemann surface of arbitrary genus g = r. The second sum is for genus case 
with more than two points specified. So we use s — r + 2 for the number of points. 
In the third sum, q = g > is the genus and r + 2 — 2q = s > 2 is the number of 
points. Thus (pTq) is equal to 



(3.17) 



s „8—2 

(*-!)(« -2) 



+ EEf-r^d - 2. 9) . „• . 

g=l s = 2 ^ ^ >' ' 

where we used Euler's formula 

C(l-2 5 ) = -^, 

and the fact that bo = 1 and &2<j+i = for q > 1. Note that the first two summations 
of (3.17) are actually the special cases of the third summation corresponding to 
s = 1 and g = 0. Thus we have established: 

Theorem 3.2. 



lim log.4 f — / exp -V ■ 



trace(X J ') dfi(X) 



y (2 g + 5 -3 ) !(2 ?) ( 2g -l ) s 2+2g+ 

2-2g-s<0 
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Since the asymptotic expansion is unique, from ( p.5[ ) we obtain 
(3.18) 

(-1)^ (2g + s-3)l(2g)(2g-l) 



E 



|Aut(r 



rib ) 



Connected ribbon graph T r n, 
with vertices of degree >3, 
X(r)=2-2g-sM^Tib)=s 

for every g > and s > subject to 2 — 2g — s < 0. 



(2g)\s\ 



C(l " 2 S ) 



Remark. If we have taken into account the values of c(n) and N in the above com- 
putation, then we will sec that all the constant terms appearing in the computation 
automatically cancel out. 

Let us examine a couple of examples. 

Example 3.1. The simplest case is g = and s = 3. The underlying graph V of a 
ribbon graph T r n,, whose topological type is S 2 minus three points, should satisfy 

(3.19) X(r) = u(r)-e(r) = 2-2s-s = -l and 

(3.20) 3v(T) < 2e(T). 



Eqn.(3.19) gives the Euler characteristic of a tri-punctured sphere, and Eqn.(3.20) 
states that every vertex of T has degree at least 3. It follows from these conditions 
that 

e(r) < 3. 



There are only three graphs in this case, as shown in Figure 3.2 





Figure 3.2. Ribbon Graphs for g = 0, s = 3 
The automorphism groups of these ribbon graphs are 62 x Z/3Z = S3, Z/2Z, 



and again Z/2Z, respectively. Thus the left hand side of ( 3.18| ) is 

(-1) 3 (-1) 3 (-1) 2 _ 1 
3! 2 2 6' 

The r ight hand side is coming from the term n 3 (— z) 1 of the second summation in 
( 3.17 ). The value is, of course, 

1 1 
6' 



3(3-l)(3-2) 



It can be also computed from ( p.lfe 



(2ff + S -3)!(2 g )(2ff-l) (2g)(2g - l)b 2g 

{2g)\s\ U g) (2g)m(2g) 
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Example 3.2. The next simple case is g = s = 1. Since the Euler characteristic 
condition is the same as in Example 3.1, the only possibilities are again graphs 
with 1 vertex and 2 edges or 2 vertices and 3 edges. There are two ribbon graphs 
satisfying the conditions: Figure 2.5 and the graph on the right in Figure [2.10 . 
The first one has 62 x Z/3Z as its automorphism group, which happens to be a 
degenerate case of the semi-direct product. The automorphism group of the second 
graph is Z/4Z, as noted in Section |[ Thus we have 

til! + (zl)! = l = _ c( _ 1) . 

6 4 12 ^ V ; 

4. KP EQUATIONS AND MATRIX INTEGRALS 

There are no analytic methods of evaluating the Hermitian matrix integral 

dfi(X) 



Z n (t,m) 



J exp ^-itrace(X 2 )^j cxp ftrace^ — X 3 ^ 



N 



However, there is an interesting fact about this integral: it satisfies the system of 
the KP equations. In this section we give a proof of this fact. 
To investigate the most general case, we define 



(4.1) 



Z n (t,m,4>) 



f ( \ 

/ exp trace > — X J 



dfi(X) 



where <j>(X) is a U (n)-invariant function on Ti n which is determined by n functions 
4>o(k), • • • , n _i (fc) in one variable in the following manner: 

(4.2) 

( (f>o{k ) 4>i(k ) ... ^n-l(fco) ^ 



det 



4>(X) =4>(k ,ki,---k„-i) = 



\<Po(k n -i) 4>l(k n -l) 



4>n-l(k n -l)/ 



A(k ,ki ■ ■ ■ ,kn-i) 



where ko, ■ ■ ■ ,k n -i are eigenvalues of X. Unlike (2.1), we allow terms containing 
t\X and tiX 2 in the integral (4.1). Using (3.6), we have 



Z„(t, m, ( 



J exp ^trace ^ <l>( x ) 



dfi{X) 



N 



f n—1 2m 



N 



I ex p EE- fc " A(fc ,... ,&„_!) 

JR n \ ._n 1 01 



det (<pj(ki)) dko ■ ■ ■ dk n - 



v i=0 a=l 



Here we need a simple formula. Let <po{k) 1 
2n arbitrary functions in k. Then 



(4.3) 



,(f> n -i(k) and tpo(k), ■ ■ ■ ,ip n -i(k) be 
det [(t>i(ke)] ■ det [ipjih)} = det [^(fc^yj) • ipj(k aU) )] , 



cre&r. 



where a runs over all permutations of & n . To prove ( |4.3| ), we calculate the left 
hand side by the usual product formula of the determinant. Then it becomes a 
summation of n n terms. Because of the multilinearity of the determinants, only n! 
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of these terms are nonzero. Rearranging the n\ terms, we obtain the above formula. 
Using this formula for ipj(k) = fc 3 , we obtain 

Z n (t,m,<f>) 
_ c(n) 



/ ex P X X ~ k i X dct (<M fc <^)) fc <^)) dk ■ ■ ■ dfc„_i 

(n) f /n—l 2m ^ \ 

/ Yl exp 5m TT^W dct )*£(<)) dk ° ■ ■ ■ dkn ~^ 

jRn it£6„ \i=0q=1 / 

X] det ( / ex p( XI ) ^'(^w) fc a(t) )*>••• 

o-ee„ \ Jr " \q=i / / 



A r 



c(n) 
TV 



c(n) 
iV 

c(n) 
AT 



det ^ y exp ^ ^ 

det^y exp^X . 



The above computation makes sense as a complex analytic function in 
(*!,••-, t 2m -i,t 2m ) G C 2 ™- 1 x {t 2m 6 C|i?e(< 2m ) < 0}, 

on which the integral converges, provided that grows slower than exp(fc 2m ). 

To compare our i^-'s with the standard time variables in the KP theory, let us set 

rp _ ta 
l a — 

a 

Now we use the formula 

(2m \ oc 

a=l 

where 



r=0 



rpTi\ rpU2 rrins rpTl^r. 

1 1 ' 1 2 ' 1 3 ' ' ' ± 2m 

nil ■ n 2 l ■ n 3 l ■ ■■n 2m \ 



(4.5) p r (T) = J2 

ni+2n 2 +3n 3 -{ h(2m)n 2m =r 

is a weighted homogeneous polynomial in Q[Ti,- • • ,T 2m ] of degree r. The rela- 



tion (4.4) holds as an entire function in Ti, • • • ,T 2m and k. Note that we have 



encountered this formula already as (l.£). From (4.4), we have 



Z n (t,m,<j>)=n\-^-det (^J ^ Pr (T)k r ^(fytfdkj 

( r°° ^ 
■det / ^pr-iiT^faifydk 

\ J — CO n , 



= n\ ■ 



N 

where we define p r (T) = for r < 0. 

Lemma 4.1. Let 4>j{k), j = 0, • • • , n — 1, fee a function defined on K smc/i £/ia£ 
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exists for all r > 0. Then as a holomorphic function defined for Refom) < 0, we 
have 



Pr -i(T) / k r ^{k)dk 



as t 



2m 



Proof. The argument is the same as the one we used in Section [y. We choose a 
fixed ti m so that Reifim) < 0- Because of the uniform convergence of the power 
series expansion of the integrand, we can interchange the integral and the infinite 
sums for a = 1, • • • , 2m — 1. Using fll.8|), (4.4) and (4.5), we have 



I-oo /2m \ \ 

J expf 5^T o Jb a j0 J -(jfc)Jfc < dfcJ 



= E 



711=0 

oo 



ni! 



00 /t~i H2r, 



"2m =0 



^ /-C 

E^( T ) / 

„ — n J— c 



k i+r (j>j{k)dk. 



r=0 



Thus we have established 



□ 



.A m, 0)) = n! 



(4.6) 



where 



c(n) 



(oo -oo \ 

^p r _i(T) / k r <f>j(k)dk 
r=o "'- 00 y 



6 



det |^p r _;(T)£ rj J , 



TV 



k r (j} 3 (k)dk. 



We recall that the determinant in (4.6) is an n x n determinant. Sato |12j proved 
that any size determinant of the form 



(4.7) 



det ^2p r -i(T)£ rj 



satisfies the Hirota bilinear form of the KP equations. He also proved that every 
power series solution of the KP system should be written as (4.7), allowing certain 
infinite determinants. A necessary background of the KP theory can be found in 



We have thus proved the following theorem. 
Theorem 4.2. If 4>j(k), j = 0, • • • , n — 1, satisfies that 

k r (j)j (k)dk < +oo 

for all r > 0, then the asymptotic expansion of the matrix integral Z n (t,m,(f>) 
satisfies the KP equations with respect to Ti,T%, ■ ■ ■ ,T2 m . Moreover, if we choose 
a value ofT^m such that Re(T2 m ) < and fix it, then Z n (t, m, <fi) itself is an entire 
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holomorphic solution to the KP equations with respect to (T\,T2,- 
C . In particular, 



,T 2 , 



l) e 



J?? 



m(Ti,T 2 ,T 3 , ■ • • ) = 55^2 log(^ n (i,m,^)) 



is a meromorphic solution to the KP equation 



;U22 



1 



"3 



-uin — 2>uu\ 



where Uj denotes the partial derivative of u with respect to Tj . 

The formula we have just established is a continuum version of the famous Hirota 
soliton solution of the KP equations |Q . The most general soliton solution of the 
KP equations due to Mikio and Yasuko Sato depends on nM + M parameters cy 
and Xi , where < i < M — 1 and < j < n — 1 . Let 

2m 

V (T,k) = J2 T <xk a . 

a=l 

Then Sato-Sato's soliton solution is given by 



0<i o <---<i„-i<Af-l 



1 ^it> 



.Jdet 



\Cin-iC 



Cz n— 1 



This coincides with our Z n (t, m, 4>) if we take 



A/-1 
i=0 



Therefore, our matrix integral Z n (t, m, 4>) of ( }4.1| ) with (4.2) is indeed a continuum 
soliton solution of the KP equations. 

So far we have dealt with the matrix integrals with a fixed integer m in this 
section. As before, we can take the limit m — > 00 of these integrals, which gives 
formal power series solutions of the whole hierarchy of the KP equations. Note 
that the determinant expression of (4.6) does not have any explicit mention on the 
integer m. Therefore, we have obtained the third asymptotic formula for the matrix 
integral: 



(4.8) 



lim A {Z n (t, m, (j>)) = n\ 




(T) / k^ji^dk 



5. Transcendental solutions of the KP equations and the 

Grassmannian 

There are several different ways to construct solutions to the KP equations. The 
Krichever construction and its generalizations are based on the correspondence 
between certain points of the Grassmannian of Sato |l2| and the algebro-geometric 
data consisting of an irreducible algebraic curve (possibly singular) and a torsion- 
free sheaf on it 0. These solutions deserve to be called algebraic, because they 
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carry geometric information of algebraic curves. Let us call a solution to the KP 
equations transcendental if no algebraic curve corresponds to this solution. The 
natural question we can ask is: how can we construct a transcendental solution? 

In this section we show that the Hermitian matrix integrals we have been dealing 
with in the earlier sections are indeed transcendental solutions. 

The technique we show that these matrix integrals are transcendental solutions 
is based on the observation that the points of the Grassmannian corresponding to 
these solutions satisfy a peculiar sl(2) stability condition. Since these solutions are 
deeply related to the moduli theory of Riemann surfaces, the appearance of sZ(2) 
is mysteriously suggestive. At present we do not have any geometric explanation 
of the relation between the KP equations, the sl(2) stability on the Grassmannian, 
and the moduli theory of pointed Riemann surfaces. 

Let V = C((z)) denote the field of formal Laurent series in one variable z. We 
fix its polarization 

(5.1) C((z)) = C[z- 1 ]0C[[z]]-z. 
For a vector subspace W C V, there is a natural map 

(5.2) -/w-W^V — > V/C[[z]]z = Clz' 1 }. 
The infinite-dimensional Grassmannian is defined by 

(5.3) Gr = {W C V | 7w : W — ► Cfc -1 ] is Fredholm of index 0}. 

The big- cell of the Grass man nian is the subset of Gr consisting of vector subspaces 
W C V such that 7w of (5.2) is an isomorphism. 

Let W be a point of the big-cell of the Grassmannian. We can choose a basis 

(w ,wi,w 2 ,- ■ ■} 

for W such that 

oo 

(5.4) Wj =z~ j + Yicijz*, 5=0,1,2,-... 



The Bosonization is a map 

(5.5) Gr— fPfq^Ta.Ta, •••]]) 

that assigns a r-function tw to each point W of the Grassmannian. For a point 
W of the big-cell with a basis ( |5.4| ), the Bosonization has an infinite determinant 
expression 

(5.6) t w = det ^pi-j (T) + ^JVfi(T)c w J . 

The infinite determinant gives a well-defined element of C[[Ti, T 2l T3, • • • ]] in the 
same manner as we have explained in the earlier sections. Sato's formula (4.7) 
gives another expression of the Bosonization map. For more detail, we refer to M 
and|. 

The commutative stabilizer of W £ Gr is defined by 

(5.7) A w = {a G C((z)) I a ■ W C W}. 

The key idea that connects the KP equations and algebraic curves is that the 
commutative stabilizer is the coordinate ring of an algebraic curve. If the greatest 
common divisor of the pole order of elements in Aw is 1, then the Bosonization tw 
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of W is essentially the Riemann theta function associated with the algebraic curve 
C whose coordinate ring is Aw |l , ■ 

Definition 5.1. A solution of the KP equations tw is said to be transcendental if 

(5.8) A w = C. 

Remark. It is known that if Aw ^ C, then the Bosonization tw of W is a solution 
to the KP equation corresponding to a vector bundle J- on an algebraic curve C 
such that 

(5.9) H C> {C 1 T)=H 1 {C : J : )={) 

|6|. Conversely, there is a solution corresponding to an arbitrary torsion-free sheaf 
T defined on an arbitrary (possibly singular) algebraic curve C satisfying (5.9). 
None of these solutions are transcendental. 

The Hermitian matrix integral we have discussed in Section ^ gives a transcen- 
dental solution to the KP equations. 

Theorem 5.2. Choose arbitrary positive integers k and n, and let 

a = (oi,o 2 ,- • • ,02fc) € <& 2k 
be a complex vector such that Re(a2k) < 0. Define a formal Laurent series 
oc / -oo /2k \ \ 

(5.10) w 3 = J2y] Ar+ ' ex p [J2 dX J z ' r+1 ~ n e C (W) 

for j = 0, 1, 2, • • • , n — 1, and Zet 

(5.11) W(a) = {wo,w u --- ,w n - 1 ,z~ n ,z- n - 1 ,---) e Gr 

be a point of the Grassmannian spanned by Wq, W%, ■ ■ ■ , w n -i, and z~ n , z - ™ -1 , 
Then the r-function corresponding to W(a) is given by the asymptotic expansion of 
a Hermitian matrix integral: 

(5.12) 

T w(a) = nm -4 / exp TjtracefX^) exp a„trace(X p ) ) dX 

where we take i?e(T 2m ) < /irsi one? t/ien Zei m — > oo to determine a well-defined 
formal power series in C[[T 1 ,T 2 ,T 3 , •■■]]. Define a linear differential operator 

(5.13) L i (a) = z 1 + - ~ z + 2^M<V* 

for i = —1, 0, 1. These differential operators satisfy the sl(2, C) relation 

[L l (a),L :j (a)} = (i - j)L i+j (a). 

The point W(a) of the Grassmannian satisfies the non- commutative stability con- 
dition 

(5.14) Li(a) ■ W{a) c W(a), i = -1, 0, 1. 
Moreover, Tw( a ) is a- transcendental solution of the KP equations. 



:!8 



MOTOHICO MULASE 



Proof. The function 

/2k \ 

exp I a A1 trace(X M ) ) 

is a special case of the function (j)(X) defined in (4.2). Thus the results of the 
previous section proves that T\y( a ) is a r-function of the KP equations corresponding 
to the point of the Grassmannian W(a). 

Let us first prove that the sl(2) stability condition (5.14) implies that the com- 
mutative stabilizer is trivial: 

Aw(a) = C. 

Suppose f(z) e Aw( a ) C C((z)), and let ord(/) = v > 0, where we define the ■pole 
order by 

ord^) = v. 
Since L_i(a) and / stabilize W(a), 

df 

[L-i(a),f] = z 2 -j- € A W (a) 
also stabilizes W(a). Note that 

ord([i_i(o),/])=i/-L 
Thus we can immediately conclude that 



A W {a) = C^ 1 } 



But then 

2k 



- ,- , v 1-/, 2 d (3n — 1) — (n — 1) 
(5.15) L-^a) -Y^n.z ^ = z j- z + - -* 

pe=l 

stabilizes W(a). Since the new stabilizer ( 5.1 5| ) decreases the order of elements of 
W(a) exactly by 1, W(a) must have an element of arbitrary negative order. But 
this contradicts to the Fredholm condition of W(a). This means Aw( a ) = C, hence 
T w(a) is a transcendental solution. 

Now all we need is to show (5.14), which can be verified by a straightforward 
computation. First, we note a simple formula 

2 k 



/oo / 2fc \ 

aA Q - 1 exp a ^J dX 

/oo / 2k 

-oo i \ i 



M A M dX. 



Let us compute the effect of the differential operators ( 5.15 ) on the basis elements 
of W(a). First, we have 



L_i(a)wj = I z 2 1- nz ■ 

\ dz 



2k \ oo />oo 
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= ]T(r + l)z r+2 ~ n / A r+ ^'eS?=i a ^dX 

r=0 

r=0M=l 00 
°° /■ 00 

= J2 rz r+l-n A r +j -l e Ef =1 <vV rfA 

2fc M-2 

M=l r=0 "'- 00 

2 OO OO 

, 1 , — 1 J— 00 



1 r—fi—l 



00 /-oo 
r=0 

2k OO /»CO 

M=lr=0 

°° /'OO 

r=0 

2/C M-2 „ QQ 

+ EE Z, ' +2 ^ / A r+ >^eS?=i^A" rfA 

°° />oo 
_£y+l-„ / (r + j)A r +J -l e E? =1 a,A^ A 

r=0 "'- 00 

2k fi—2 „ qq 

= -jwj-i + z r+2 ~ n -» / A r+ >a^eSf=! a ^ M rfA 

1 n J —OO 



f_i=l r =0 

E W(a) 



for all j = 0, 1, 2, • • • , n — 1. Note that u>_i does not appear in the above compu- 
tation because of the combination jwj-i. For the basis elements z~ n , z~ n ^ 1 , ■ ■ ■ , 
we have 



L_i(o)2- n -'= lz 2 ^-+nz 



dz 

2k 



2k \ 
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for all i > 0. We note that the term z~ n+1 does not appear in this computation. 
Thus we conclude 

L-i(a) • W(a) C W(a). 

For j = 0, we have 



/ d 3n - 1 



2 



2fc \ oo ->oo 

M=l / r=0 "'- 00 



jr (r + Il + V + i-n / A^e^-^dA 

oo 2fc /-oo 
n i J —oo 



r=0 fi=l 

n + 1 



1 , , roo 



o 

2fc (U— 1 



11=1 r=0 

2 /c oo „ qq 

££y + i- n -„ / A^V^e^^^^^A 

. n J —oo 



2 

2fc (U— 1 



2fc OO /■ oo 

H=lr=0 J- 00 

+ ^rz r+1 - n X r+ ie^l=i a ^dX 

„ — n J —oo 



1 00 

n + 1 



2 J 

2k n-1 . 



^ poo 

+ EE zr+1_ "^ / y +j ^e^ a ^dX 
~ S zr+1_ " / ( r + - ? + 1 )^ r+i e E ''= 1 a " X "dX 

r=0 

' H=lr=0 



x "dX 



€ W(a) 

for all j = 0, 1, 2, • • • , n — 1. It is obvious that 

Lo(o) • z- n -* g W(a) 
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for i > 0. Finally, for j = 1, we have 
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Li(a)Wj 



(, 2fc \ oo ^oo 

+ (2n - IK 1 + ^ ^r"" 1 53 z r+1 -" / A r+ ^?=i a " AM rfA 

°° /'OO 

= 53 + »)« r ~" / A r+J 'e E '=i a * A "dA 

oo 2 Aj /* qq 

+5353^-"-" / v+v^^^^dA 

, — n ., — 1 J —oo 



r=0 n=l 
oo /.oo 



= ^ (r + n + l)z r+1 -" /°° A r+ '' +1 e£?=i»M v * dA 

r=-l "'- 00 
2/c „qo 

+ E E ^"""^ / A'+V^" 1 *^d\ 



2k oo . v 

+ 53 53 z^""" / X r +^a^l^ a ^dX 

J — rm 



/ — OO 
/■OO 



,— 1 r— 



/oo °° ^00 

AJ e l£=i <v A "dA + 53 (r + n + l)z r+1 -™ / A r+J+1 e^=i a * A "dA 
-oo r= Q J —oo 



2k n 



-53 ^V"""" /°° A^V^-^A 

p=lr=0 



2fc oo 



/OO 00 <-00 

AJ'eS^i <v A "dA + 53 (r + n + l)z r+1 - n / A r+J+1 e^=i a " x " dX 
-OO r _Q J —oo 

2k fJ> „qo 

+ E E * r " n "" / A'+V^" 1 ^ A "dA 



f_i=l r=0 

/>oo 



- ^ z r+1 - n / (r + j + 2)\ r +i +1 eZT=i a ^d\ 

r=0 "'- 00 

/oo 
A^^i a ^dA + (n - j - lK+i 
- OO 

2/c M /.oo 

+ EE zr_ " _AI / A^+Va^'i^dA 

(i=lr=0 

G W(a) 
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for all j = 0, 1, 2, • • • , n — 1. Note that the term w n does not appear in the compu- 
tation. It is again obvious that 

L x (a)-z- n - 1 G W(a) 

for i > 0. This completes the proof of the sl(2) stability of W(a), and hence we 
have established the theorem. □ 

The action of these sl(2) generators on W(a) is very subtle, and it does not seem 
to allow any generalization. For example, the above proof does not apply for the 
Virasoro generators Li(a) other than % = — 1,0, 1, although the operators L,(a) are 
defined for alH e Z and they satisfy the Witt algebra relation 

[Li(a),Lj(a)] = (i-j)L i+J (a) 

for i, j e Z. 
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